Heavily Enhanced Dynamic Stark Shift in a System of Bose Einstein Condensation of 
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The dynamic Stark shift of a high-lying atom in a system of Bose Einstein condensation (BEC) 
of photons is discussed within the framework of nonrelativistic quantum electrodynamics (QED) 
theory. It is found that the Stark shift of an atom in BEC of photons is modified by a temperature 
dependent factor JbeCj compared to that in a normal two-dimensional photonic fluid. In photonic 
BEC, the value of Stark shift is always greater than that in two-dimensional free space. Physical 
origin of the phenomenon is presented and potential application is also discussed. 



o 



<Z3 

^— > ' 

c : 

=5 
+-» ! 



i 

C 

O 
o 



> 
m 

o 



X 



I. INTRODUCTION 

Bose Einstein condensation (BEC). predicted by Ein- 
stein [l| in 1924, is a remarkable state of matter. For 
a bosonic system at low temperature, significant friction 
of particles would spontaneously condensates into zero 
momentum state with minimum free energy of the sys- 
tem. BEC was first observed in rubidium by Anderson 
et al. [2J in 1995. Thereafter, the phenomenon was re- 
vealed in many other systems of atoms or quasi-particles, 
such as sodium Q , lithium [1] , cesium Q , potassium @ , 
hydrogen polariton Q, etc. However, BEC of the 
simplest boson system, the photon system, was not ob- 
served until 2010 What impeded us so long from 
realizing Bose Einstein condensate (BEc, lower case "c" 
to distinguish from BEC) of photons is that in a nor- 
mal blackbody radiation cavity, the photon is massless 
and the photon number is not conserved which leads to 
a vanishing chemical potential. In 2010, Klaers et al. @ 
have overcome both obstacles by confining laser pump- 
ing light in a two-dimensional microcavity which is filled 
with dye and bounded by two highly reflective concave 
mirrors. They established the conditions required for the 
light to thermally equilibrate as a gas of conserved par- 
ticles rather than as an ordinary blackbody radiation. 
As is well known, quantum optical effects of atoms are 
not only dependent on their internal structure, but also 
on external electromagnetic environment. People have 
explored various systems in modified electromagnetic en- 
vironment such as in dielectric medium [Tol. , photonic 
crvstals(T2l [l3|. and optical microwave guides 14]. The 
results show that the energy-level shift is modified ac- 
cordingly. For example, Wang et aZ.[l5[ predicted that 
the dominant contribution to the Lamb shift comes from 
emission of real photon in photonic crystals that the 
Lamb shift can be enhanced by 1 or 2 orders of mag- 
nitude, termed as 'giant' Lamb shift. In a recent work, 
we investigated the Lamb shift of a hydrogen atom in- 
side a Kerr nonlinear body (KNB) and also found that 
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the modification of the Lamb shift is a 'giant' one [la ]. 
What is more, we found the dynamic Stark shift in KNB 
revealed a temperature dependency fl~7j . Inspired by the 
experimental demonstration of BEC of photons, in this 
paper, we aim to investigate the dynamic Stark shift of 
an atom in a system of BEC of photons. We find that 
in such a system, the dynamic Stark shift of a high-lying 
atom can be heavily enhanced. 

The remainder of this paper is organized as follows. In 
Sec. II, we configure the microcavity for BEc of photon 
system and theoretically model the BEC of weakly inter- 
acting photons in a two-dimensional optical cavity. The 
dynamic Stark shift of a high-lying atom in photonic BEc 
is derived in Sec. III. In Sec. IV, we discuss the order 
of modified shift and the role of dimension. Finally, we 
make a brief conclusion in Sec. V. 



II. BOSE EINSTEIN CONDENSATION OF 
PHOTON SYSTEM 

In this section, we theoretically establish the model for 
BEC of a two-dimensional photon gas. We first base the 
design on Klaers' work to configure the two-dimensional 
microcavity in which photonic BEC may realize. Then 
we turn to reexamine the dispersion relations of two- 
dimensional free and weakly-interacting photon gas. Fi- 
nally, we discuss the conditions of BEC of photons. 



A. Microcavity resonator 

The first difficulty stands in our way is that the chem- 
ical potential of photon system is vanishing in normal 
blackbody radiation, that is, the photon number is not 
conserved when the temperature of the photon gas is var- 
ied. However, a prerequisite of BEC is that the chemical 
potential be non- vanishing [18j . For a black body radi- 
ation, photons disappear in the cavity walls instead of 
occupying the cavity ground state Q. Luckily, number 
conserving thermalization of two-dimensional photon gas 
was experimentally achieved by Klaers et al. @, theo- 
retical work was previously done by Chiao et al. 
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Typically, a microcavity resonator is required to be corn- 




High reflectivity mirrors 



FIG. 1. Schematic of the optical microcavity. Consisting of 
two high-refiectivity curved mirrors which are fulfilled with 
dye solution. An atom is immersed inside the microcavity 
away from the laser path. 

posed of highly reflective mirrors with reflectivity over 
99.997% at both end, acting as the walls, as shown in 
FIG. 1. Separation between surfaces of adopted curved 
mirrors are described by equation D (r) = Dq — 2(i? — 
V R 2 — r 2 ) with r being the distance from optical axis. 
A laser of 532 nm pumps dye onto the cavity at an an- 
gel of 45° to optical axis. The atom should be located 
inside of the cavity and completely immersed in solution 
but away from pumping beam to avoid possible pumping 
light excitation. Generally, the atoms couples with pho- 
tons forming polaritons, a bosonic quasi-particle. How- 
ever, due to the frequent collisions of photons and dye 
solution, this coupling is broken pol [2l|. Thus we may 
convince ourselves that the BEC we observed here is real 
BEc of photons instead of that of polaritons. 

B. Free Photon Dispersion Relation 

The mirrors have confined the photons to a two- 
dimensional plane, simultaneously imposing a boundary 
condition on both sides of inside electric field. The al- 
lowed z-component which is along the optical axis of wave 
vector of photon is k z (r) = n'n / D(r). n' is thus an in- 
teger. The electric field vanishes at the reflecting surface 
thus imposing a quantization condition for a component 
of photon wave vector which denotes as z. Photons are 
trapped inside the resonator. In this two-dimensional 
plane, a plane-mode wave propagates at small angle with 
respect to z axis, the effective mass of a non-relativistic 
particle is m 7 = hui/c 2 22]. 

For a two-dimensional noninteracting photon inside the 
cavity resonator, the de Broglie dispersion relation or 
energy-momentum relation under small-angel approxi- 
mation is 

p 2 

e(p) = m y c 2 + ——. (I) 



We note that in present cavity resonator k z , which is 
greater than that of x or y component k x , k y under the 
small angle propagation. Thus, we ignore the contribu- 
tion of k z , giving 

2 2,2 
P = P X +Py 

Confining our sight to the two-dimensional plane, we may 
view free photons as normal particles in a transverse res- 
onant plane. This may owe to the fact that number con- 
servation is achieved by scattering off process of photons 
and dye molecules, which make it a canonical ensemble. 



C. Bogoliubov Dispersion Relation for Photon Gas 

There exists Bose Einstein condensation in an ideal 
Bose system at absolute zero temperature. This feature 
survives in the case of the weakly interacting Bose sys- 
tem, since as the interaction vanishes, one should recover 
the Bose Einstein condensate state[23|. From the work 
by Klaers et al. [9J] , we may find the magnitude of the in- 
teraction energy of photon knowing the intensity of light 
I(r) and relative refraction index /j, 

E- mt = m 7 c 2 nl(r). (2) 

We rewrite it using what applied in Gross— Pitaevskii 
method, introducing a dimensionless parameter [24j g = 
—m^c & /2'Kh?nq to describe the strength of interaction. 
Here, c is the speed of light, n is the the quantum con- 
centration of current system, and q is a constant. So we 
rewrite the interaction energy E- ln t in term of g 

E in t = J^/m^gNon 2 . 

In this cavity, g ~ (7±3) x 1CP 4 , far below that reported 
in two-dimensional quantum gas experiment, in order of 
—2 ~ — 1. This indicates that the photons are weakly in- 
teracting which enable us to apply Bogoliubov transform 
in such a two-dimensional quantum gas system. 
Since this system is open, i.e., it is connected to an ex- 
ternal reservoir of particles, so that the average of to- 
tal particle number fluctuate around an average value, 
then the total particle number needs only be conserved 
on average. To formally illustrate it in the expression of 
Hamiltonian, we introduce Lagrange multiplier (iN open , 
like that done in statistical mechanics. 

H = H' — M iV open (3) 

where iV p Cn = ^ p o-\,o,p is total number operator and 
fj, the chemical potential. Experimentally, pump beam 
in the configuration is the physical implementation of 
external reservoir. The overall energy is thus given by 
the equations above. Now, we study the Hamiltonian of 
the system in details 

H = H int + -ffnon-int- (4) 
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For a free boson system, with de Broglie dispersion rela- 
tion in Eq. ([T]), the Hamiltonian is 



-ffnon-int = ^ e(p)a p a p . 



(5) 



We introduce the creation and annihilation operators ap 
and a p for photons with momentum p, respectively. We 
know that in a weakly-interacting Bose gas at absolute 
zero temperature, most particles with zero momentum 
occupying condensate state, are in great number N . For 
such ground state \Nq), the zero-momentum operator 
and cio obey the following relations 

a \N ) = 



where 



and 



e'(0) = JV e(0) - ^V N 2 



e'(p) = c(p) + N V P + N Vo - (i 



(10) 



(11) 



are modified photon energy. 

Following Bogoliubov[25j, we now diagonalize the Hamil- 
tonian, that is, transforming the photon operators into 
quasi-particle ones. Since the Bose commutation relation 
preserve in such case, transform can be introduced and 
applied in the following form 



4 \N ) = VN + 1 |JV ) • 

Noting that Nq is large enough in this case, so, we may 
approximately take ^/Nq + 1 « \/Nq. The operators 
also satisfy Bose commutation relations 



where 



b p = Updp + v p a_ p , 
b p = Updp + v p a- p 



(12) 



[dp, a\] 



[a p ,a q ] = [a P ,at] 



0. 



(6) 



The interaction term describes the momentum transfer 
of two photons arising from potential energy V(p), it 
represents annihilation of two photons with momentum 
p and q, with the creation of two photons with momen- 
tum p + n and q — k. Specifically, the annihilation 
of two zero momentum photons giving rise to non-zero 
ones with momentum transfer k within this pair is 
the phase transition from condensate state to normal 
state. The interaction is reasonably in the following form 



= J2 \ V ( K ) a P+K a l- 



q-K, a p a q 



(7) 



In BEC, since only small momenta are concerned, we 
consider the zero-momentum term of interaction poten- 
tial 



V = 



Airft 2 



(8) 



where s — g/2Tt is s-wave scattering length and m 7 
is effective mass of photons [18j. And we assume that 
V(p) — V{— p) is valid. 

Since in BEC, particles are macroscopically occupying 
the ground state, N is relatively large, that is N w 
JVo > 1. Terms with a ao = No dominates. Under 
weakly interaction approximation, for interaction terms 
in only that with 7V and iVp are preserved, so the 
Hamiltonian H is 



H~ e'(0) + ^e'(p)a p a_ p 

p#0 



+ N Vp{apa [ _p + a p a_ 

p#0 



(9) 



and 



[b p ,bl] = S P: q, 
[b P ,b g ]=[blbl] 



0. 



(13) 



A diagonal form of weak interaction term of p and — p 
momenta photons are expected to be 



En 



(14) 



where E is ground-state energy. Eq. (fl"4"|) shows that 
the original system can be reexpressed in term of new 
set of operators b p and b^ q which respectively represents 
the annihilation and creation of a quasi-particle with en- 
ergy t{p). Their motion satisfy the following equation in 
Heisenberg picture 



ih^- = [bp,H] = ~tip)b f 



(15) 



Substituting Eqs. (fT2|) and (fT3| into original Hamilto- 
nian (Eq. ©) and comparing with Eq. (IT41 . we get the 
following diagonalization conditions 



l p v p 



NoVp/tip), 



Up = -[t'ip)/~tip) + l], 

vl = \[t'ip)/tip)-l]. 
Solving Eq. (jTTJ) and Eq. (TTSl) for tip), we have 

lipf = t\pf - N%V£ 

= tip) 2 + 2tip)N n V , 



(16) 



(17) 
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Ip 2 n v 



P' 

4m2 ' 



(18) 



This is the Bogoliubov dispersion relation for the present 
system. Specifically, at critical temperature, the two 
terms in Eq. (1181) are equal, so that critical momentum 
p c = 2y / m 1 NoV Pc [23]. At low temperature, we adopt 
long- wavelength approximation so that 



In present configuration, nq — 1.00828 x 10 15 and n — 
7.76374 x 10 18 . Condition is satisfied. 
Here, the dye solution serves as the heat bath and equi- 
librates the transverse modal degree of freedom of pho- 
tons and dye molecules via absorption and re-emission 
processes. At room temperature, photon frequency is 
above the low-frequency cut-off, in contrast to the case 
of a blackbody radiator, for which the photon number 
is determined by Stefan-Boltzmann law thus reveals a 
temperature dependency. 



= c\p\ 



(19) 



where c = y/ No Vb/m 7 is the sound velocity [T^. Actu- 
ally, the fluctuation of photon gas described by c and 
bp bp suggest that the excitations here is phonon which 
can be reexpressed as 



III. DYNAMIC STARK SHIFT 

Now, we turn to the electromagnetic environment of 
the atom inside the cavity. First, we write down the 
Hamiltonian of the system 




y/N-aT 



where N is the total number of photons in cavity and 
a = m-fkB^p^oP 2 i s a constant independent of p 
nor temperature T, but concerns the distribution of 
momenta [18| |. 



D. Condensation Condition 

We now examine the condition of BEC. All phase tran- 
sition has a critical point, here is the critical temperature 
T c . Only at low temperature, particles are occupying 
ground state. Above T c , the condensate phase vanish, 
photons are free. For a normal three-dimensional BEc, 
the critical temperature is [26[ 



3.31fi 2 N ' 



(20) 



where V is the volume of the resonator. However, in two- 
dimensional harmonic tarp, a revising factor ~ 0.78 

is applied with hard core calculations, see Eq. (18) of 
Ref. [13] ■ For photons, the degeneracy g — 2. Calcu- 
lation shows that in present configuration, the critical 
temperature is about 578.062 K (for data see the cap- 
tion of FIG. 2). That means experiment can be achieved 
even in room-temperature. However, there exists another 
condition that quantum concentration constrains 



n > n Q 



m 7 fc B r .3 /2 

1 2nh 2 ' 



(21) 



where n is the quantum concentration of current system, 
we find it using 



N ( 



27T?l 2 , 

m^k^T 1 



-3/2 



(22) 



H + H! 



(23) 



where the free Hamiltonian Hf rcc is that of a bare atom 
and H[ nt is the interaction Hamiltonian between atom 
and electric field. This atom can either be at high-lying 
(Rydberg) state or low-lying state. Typically, we con- 
sider an atom with only one electron on its outermost 
electron shell. Atoms like hydrogen and alkali like ru- 
bidium, potassium are widely used in experiment. Much 
work has already been done by Hollberg et al. and 
Zimmerman et al. (29| both theoretically and experimen- 
tally. For such an atom, we express the Hamiltonian in 
form of 



EL In) 



(24) 



The other term H has been derived in the previous sec- 
tion, see Eq. (fl4|) . Our interest now is focused on the 



H[ ai term, the interaction Hamiltonian between atom 
and photon fluid is 



E 



(25) 



where r is the radius vector of the electron, E is the 
electric field induced by BEc of photons, and e is the ab- 
solute value of the electron charge. To get the expression 
of H' int , we need the expression of E. Generally a vector 
potential of electromagnetic field is 



2Ve Q ujp 



(a p (t) + a p (t))e p . (26) 



It is worthy of noting that to emphasize the time related 
term e~ luJpt , we add "(i)" to indicate that operators a p 
is actually time-dependent. For b p , it is e _i " pt term. ui p 
is the angular frequency of phonon with momentum p. 
We have shown that the interaction between the field in- 
duced by weakly interacting photons in BEc and atoms 
can be viewed as the collection of non-interacting exci- 
tations that perturb the atoms. Using Bogoliubov trans- 
form (Eq. (H2)), we rearrange the photon BEc induced 
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electric field into excitation operators described form. 
Using equation (fT2"j) and (fT5|) . the vector potential is 



A = E 



^(6 p (t) + 6t(t))e p (27) 



where 




At critical momentum p c , the interaction potential has a 
harmonic form 23] 

V p = \p 2 (28) 

where A = l/4m 7 iVo is a constant with a unit of mass -1 . 
A reasonable assumption is adopted that A is indepen- 
dent of p under current near zero temperature approxi- 
mation. With this we evaluate the value of e' 3 



=2/3 



N V P 



N Xp 2 
p 2 /2m~ l 



(29) 



= y/l + 2m 1 N X 
= \/S/2. 

Now, we turn to the electric field E. By using Maxwell 
equation E 



-7^1 the electric field is 



E 



E 



2Ve uj p 



(6p(t) - 6+(t)) e p . (30) 



So, we have interaction Hamiltonian 



2Ve u} p 



ojp(b p {t) - & P (i))r ■ e p . 



Now, we consider the perturbation term who bares its 
first order term being zero due to symmetry of electrical 
dipole. So we turn to the second order 



AE(4>) 



ct>^\H' mt \I,^ 



E 

n' ,7,p 

\{^;\H' int \i,^)\ 2 Np 

^E* + N p e> p -(Ei + (N p ~iy p ) 

\{^^\HL t \i,^ +1 )\ 2 



= E( 



E 



E<t> + N p e' p -{Ei + (N p + l)e' p ) 
2 2Ve ujp p E$ — E 1 ± e p p 



(32) 



l(^|r-e p |J)| 



E<P -E 1 - e 



AS 



AC 



AS 



Lamb 



It is worthy of noting that in result Eq. (|32|) , two terms 
are presented. The first one is the dynamic (or AC) Stark 
shift induced by real photons, and respectively the second 
one is the Lamb shift induced by virtualphotons or , 
say, by vacuum. Mass renormalization |30| should be 
preformed to each term followed by integral ranging from 
to U. Since we are considering a non-relativistic case, 
an upper limit of energy, here U, should be applied that 
is large enough with respect to the energy difference E^ — 
E 1 [3l|]. Thus we have 

AE AC = J%^<*f.)* y P / dt \r. A 2 x 
V 2 47re 5 v he ' ^ J 1 ' 1 e* - 1 



I 



1 



r) 



-(E^-E^/knT + t (Ef- Ei)/kvT-t' 

(33) 

where t = pc/k B T, |r 0j/ | 2 = |(0|r|/)| 2 , K = J fc 2 (r)dr 
and 



AE 



Lamb 



3 e 2 Kc ,k B T 
2 Airsnc 



(34) 



\ r <t>J\ 



{E$ - Ei)/k B T - t 



"P" in the above equations denote that a principle value 
integral should be preformed instead of the normal one 
since it is a improper integral. 

By using the normal vector potential (Eq. (f2"o]) ) and its 
corresponding electric field, we can easily get the dynamic 
Stark shift in a normal two-dimensional photon gas 



A P AC 
^•"Normal 



47T£f 



r <t>,I\ 



1 



I 



(Et - Ei)/k B T + t {Et- Ei)/k B T -t' 

Comparing Eq. Q33p with Eq. (|35p , we can rearrange the 
shift in the form of 



where 



AE AC = F B EC ■ A^Normal 



3 c 



(36) 



(37) 



is the contribution factor of BEc. In this equation, we 
may read and find what BEc-related shift is different 
from a normal one. The factor -Fbec described the 
modification with respect to the Stark shift in a normal 
photon fluid A££C mal . 



IV. RESULT AND DISCUSSIONS 



From Eq. (|36|) , we can see that the Stark shift in a sys- 
tem of BEC of photons is modified by a factor Fbec com- 
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pared to that in a normal two-dimensional fluid. To bet- 
ter illustrate the physical interpretation, we can rewrite 
it in the following form 



"BEC 



3m, 



2V y/N -aT 



In order to give a numerical impression of -FbeCj we take 
the experimental data from Klaers et al. [9(. The ef- 
fective mass m 7 « 6.7 x 10~ 36 kg, the zero momentum 
interaction potential Vq « 3.2 x 1CP 33 J, the speed of 
light c w 3 x 10 8 m/s, and the total photon number TV 
no less than 77000. As a result, Fbec is a number of 
order 4^5. FIG. 2 shows the function of -Fbec to aT. 



,x10 




FIG. 2. Modification factor Fbec as a function of modi- 
fied temperature aT. Data are based on the experiment per- 
formed by Klaers et al.®. V = 47rft 2 s/m 7 « 3.2 x 10" 33 J, 
cut-off frequency w cu t-ofl ~ 3 x 10 15 Hz, effective mass of pho- 
ton m 1 ~ 6.7 x 10 -36 kg, dimensionless parameter for interac- 
tion strength g ~ 10 -3 , and the volume V 



Near absolute zero temperature, it does not vanish but 
approximate 6.6315 x 10 4 . In spite that the normal shift 
is small, the total effect is significant. So, as long as the 
fluid is in BEC state, the shift is far greater than the 
normal one. 

Up to now, one may raise the question why the dynamic 
Stark shift is so large in BEC of photons. This can be 
answered as follows. As is known, the dynamic Stark 
shift is caused by perturbation of electric field, here, con- 
tributed by real photons. The key physical function con- 
cerning the atomic QED is the density of state (DOS) of 
the photonic system. In a non-absorbing linear medium, 
the DOS is given by 



(38) 



However, in present environment, the photons are re- 
placed by quasi-particles. The DOS has to accordingly 
change to a quasi-particle described form 



We can see that DOS in BEC of photons is a function of 
the quasi-particle velocity c, thus a temperature depen- 
dent value. At near zero temperature, c is much smaller 
than c, then the DOS now is relatively much larger than 
normal one. Much concentrated quasi-particles results in 
a heavily enhanced Stark shift. 

Additionally, we consider the role of dimension of the 
electromagnetic environment on the dynamic Stark shift. 
Unlike other photon systems, a two-dimensional system 
applied here, current resonator have a non- vanishing 
chemical potential, significantly different from three- 
dimensional photon system typically encountered in 
Planck problem. This may be accounted by the repulsive 
pairwise interactions between photons in a Bose Einstein 
condensate. To find the chemical potential fi, use 



M = 



dE 
dN 



(40) 



where Eq is ground-state energy (tp \H\ tjj). Since it is in 
superfluidity, which is necessarily Bose-Einstein conden- 
sate state, the zero-momentum photons or pairs domi- 
nates, thus the other terms are dropped out. Previously, 
we have eo = hN 2 Vo that makes 



M = N V 



(41) 



which clearly shows that a non- vanishing chemical po- 
tential is actual in current resonator. The existence of 
low-dimensional BEC has already been discussed. Back 
to 1997, Mullin theoretically showed that the Bose gas 
has a phase transition at some critical temperature T c in 
a two-dimensional harmonic trap in the thermodynamic 
limit [27} . Further, he pointed out that for dimension no 
less than 1, T c follows the same expression, however, in 
one-dimensional case, there is no condensation. 
The dynamic Stark shift in a three-dimensional normal 
fluid given by Farley et al. [32] 



67T 2 £o he 



1 



1 



{E+ - Ei)/k B T - t 



Ei)/k B T + t 
)• 



(42) 



p(Q p ) = Vl>1/tt 2 c 3 



(39) 



Compare above equation with that for a two-dimensional 
case Eq. ([55]) . we may conclude that for Stark shift in 
system of different dimensions, the dimensionality con- 
cerns the exponential of temperature term k B T '/he and t 
in integral term. A higher order of dimensions makes the 
shift far greater in a relatively high temperature. How- 
ever, in a system of BEC of photons, it can be neglected. 
In rencent years, the technology of atom manipulation 
by laser has been developing rapidly. Stark shift is con- 
sidered to be suitable to compensate the Doppler shift 
of traveling atom. In a normal electromagnetic envi- 
ronment, the Stark shift is too small to be qualified. 
However, in the BEC of photons, the heavily enhanced 
Stark shift seems promising. Due to its large value, a 
wide range of frequency in Doppler compensation might 
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be possible. What is more, the heavily enhanced Stark 
shift might be used for tunable far-infrared photodctcc- 

tors nm. 

V. CONCLUSION 

In summary, we have investigated the dynamic Stark 
shift of a high-lying atom in a system of BEC of photons 
within the framework of nonrelativistic QED theory. In 
this configuration, the effective mass and the chemical 
potential for weakly interacting photons inside the mi- 
crocavity are nonvanishing. Via Bogoliubov transform, 
sound waves exists for long-wavelength disturbances of 
the system. It is found that compared to that in two- 
dimensional free space, Stark shift of the atom in BEC 
of photons is modified by a factor -Fbec • It is a monoton- 
ically increasing function of temperature T and depends 



on a few other parameters of the system. Below the crit- 
ical temperature T c , the value of Stark shift is always 
greater than that in a normal two-dimensional photonic 
fluid to an order of 4 ~ 5. Physical origin of this phe- 
nomenon is discussed and future applications is proposed 
including Doppler shift compensation and tunable far- 
infrared photodetectors. It is hoped that the predicted 
properties will be verified in physics laboratories for the 
not too distant future. 
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